Cooling of a single atom in an optical trap inside a resonator 
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We present detailed discussions of cooling and trapping 
mechanisms for an atom in an optical trap inside an optical 
cavity, as relevant to recent experiments. The interference 
pattern of cavity QED and trapping fields in space makes the 
trapping wells in principle distinguishable from one another. 
This adds considerable flexibility to creating effective trapping 
and cooling conditions and to detection possibilities. Friction 
and diffusion coefficients are calculated in and beyond the low 
excitation limit and full 3-D simulations of the quasiclassical 
motion of a Cs atom are performed. 
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I. INTRODUCTION 

A recent experiment succeeded in trapping a sin- 
gle atom with single photons inside an optical cavity and 
in monitoring the atomic motion with the resolution ap- 
proaching the standard quantum limit for position mea- 
surements. Yet a second experiment fa has likewise re- 
ported single-atom trapping at the few-photon level, al- 
though in this case the trapping potential and diffusion 
are in fact well approximated by a free-space semiclassi- 
cal theory jj) . 

One future objective for such experiments is to use 
atoms trapped in cavities for quantum communication 
purposes, with atoms serving as quantum memories and 
photons as the transporters of quantum information [Q,^) . 
While the single-photon trapping experiments provide a 
new paradigm for quantum measurement and control, 
they are nevertheless not entirely suitable for the pur- 
pose of distributed quantum networks where qubits will 
be communicated among quantum nodes. The reason 
is the short trapping life time of the atoms as well as 
limited operation flexibility A better strategy might be 
to use the cavity QED field for quantum state entan- 
glement and distribution while an additional (external) 
trapping mechanism provides the necessary confinement 
of the atomic center-of-mass motion. For instance, in an- 
other recent experiment from the Caltech group || , mean 
trapping times of ~ 28ms (as compared to mean trapping 
times of < 1msec in the experiments ^J,^|) were achieved 
by employing a far-off resonant trapping (FORT) beam 
along the cavity axis. In that experiment the trapping 
lifetime was limited due to intensity fluctuations of the 
intracavity FORT beam 0. Here we consider the situa- 
tion of current improved experiments B in which a single 



atom is held inside an optical cavity in a stable FORT 
beam of minimum intensity fluctuations. 

Several mechanisms for cooling inside optical res- 
onators have been discussed before p-pd|. Here we dis- 
cuss in detail how the combination of an external trap- 
ping potential and the cavity QED field adds flexibility 
in predetermining where and to what degree atoms will 
be trapped and cooled. Moreover, our calculations go 
beyond the weak driving limit discussed in JhJ. That 
is, we allow the "probe" field driving the cavity to be so 
strong as to appreciably modify the dynamical behavior 
of, rather than merely probe, the atom-cavity system. 

This paper is organized as follows. In Section II we 
describe the physical situation of an atom trapped in an 
optical potential and strongly interacting with a cavity 
QED field. We give the evolution equations for both in- 
ternal and external atomic degrees of freedom and for 
the quantized cavity mode. Section III contains an ex- 
position on how we calculated friction and diffusion co- 
efficients from the forces acting on the atom. Section IV 
contains the main results of this paper: we discuss sim- 
ple pictures for cooling mechanisms, based on the dressed 
state structure of the atom-cavity system, and give nu- 
merical results for the typical cooling and diffusion rates, 
and hence "temperatures" for single atoms under vari- 
ous trapping conditions. We also study the saturation 
behavior under strong driving conditions and perform 
simulations of the full 3-D motion of atoms trapped in 
particular wells that show how the probe field transmis- 
sion is correlated with the atomic motion and how trap- 
ping times can be prolonged by strong cooling. Section 
IV F concludes with a brief discussion of a slightly differ- 



ent trapping scheme. The summary highlights the main 
results. 



II. DESCRIPTION OF PROBLEM 

We consider a single two-level atom coupled to a single 
quantized cavity mode and coupled to a (classical) far-off 
resonant trapping (FORT) beam. In most of the paper 
we assume that the FORT shifts the atomic excited state 
|e) up and the ground state \g) down by an amount Sf(t) 
(i.e., the energy of the ground state is E 9 — Sf, that of the 
excited state E e + Sp), as this is the situation purs ued in 



previous and current experiments [|6||8j]. In Section IV F , 
however, we will also study the different situation where 
both ground and excited states are shifted down by Sf 
(see e.g. |U2|). The FORT beam coincides with one of the 
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longitudinal modes of the cavity and its wavelength Xf 
is longer than that of the main cavity mode of interest 
for cavity QED, Ao- In fact, in the experiments the 
cavity length L is 104A /2 = 102A F /2. 

The position-dependent AC-Stark shift due to the 
FORT field is of the form 

S F (r) = So sin 2 (fc F z) exp(-2 / 9 2 /wg), (1) 

with S*o > the maximum shift, kp — 2it/Xf the wave 
vector of the FORT field, wq the size of the Gaussian 
mode of the cavity, while z and p give the coordinate 
along, and the distance perpendicular to, the cavity axis, 
respectively. The quantized cavity mode is assumed to 
have the same transverse dimensions woQ so that the 
atom-cavity coupling is determined by 

9(r) = 9osm(kz)exp(-p 2 /wl), (2) 

with go the maximum coupling rate and k = 2tt/Xq the 
wave vector of the cavity mode. Under conditions where 
the cavity is not driven too strongly, the atom will be 
trapped around the anti-nodes of the red-detuned FORT 
field. Thanks to the fact that Ao ^ Xf, the atom will ex- 
perience a different coupling strength to the cavity mode 
in each different well. Figure [l] shows the axial pattern 
arising from the FORT and cavity fields. For illustrative 
purposes we choose here (and in the rest of this paper) a 
cavity of length L — 16Ao = 15Af- This does not influ- 
ence the basic physics involved: in particular we note that 
the precise value of Xf is largely irrelevant on the time 
scales considered here, as the FORT field is detuned far 
from atomic resonance. The choice of L — 16Ao — !5Xp 
just means that only 8 wells out of 30 are qualitatively 
and quantitatively different. 
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1 It is in fact the Rayleigh ran ges of th e beams that are iden- 
tical, so that w% ORT /w^ v = y/\o/X F « 0.99 



FIG. 1. The FORT-induced shift of the ground state on 
axis (p = 0) in the case where 5*0/(271) = 10MHz and the 
cavity-QED coupling rate g (dashed curve), with go = 35o, as 
functions of position along the cavity axis measured in units 
of the FORT wavelength Xf, with z = at the left cavity 
boundary. The cavity length is L = 15Af- 

This is illustrated in Fig. || where we plot the value 
of the cavity QED coupling g at the anti-nodes of the 
FORT (i.e., the bottom of the trapping potential). In 
particular, there are 2 anti-nodes in which g = 0, and 4 
in which \g\ attains it maximum. 
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FIG. 2. The values of \g\/go at the locations of the 
anti-nodes of the FORT, i.e. at the points z n — (n/2 — 1/4)Af 
for n — 1 . . . 30. There are 8 quantitatively different wells. 

The cavity is driven by an external classical field 
£{t) = £q exp(wjpt), at a frequency uj p , which is used 
to probe the atom-cavity system and which may cool the 
atom at the same time. In the following, the strength 
of the driving field is indicated by the number of cavity 
photons iVe that would be present if there were no atom 
in the cavity, rather than by E§. This closely follows 
the experimental procedure for determining the driving 
strength. The relation between the two is 



with A c = lj c — Up the detuning of the probe from the 
cavity frequency uj c = kc. The Hamiltonian for the inter- 
nal atomic degrees of freedom and the quantized cavity 
mode is, in a frame rotating at the probe frequency u> p , 
given by 

H = hA c a + a + hA a a + <7~ + 2hS F {r){<7 + o~ - 1/2) 
+M (a + + a) + hg(r)(a+<7- + a+a). (4) 

Here A a = oj a — u> p is the detuning of the atomic reso- 
nance from the probe frequency. In all numerical exam- 
ples given below the cavity frequency is chosen to coincide 



2 



with the atomic frequency, so that A c = A a . The quan- 
tity A p = — A a is then referred to as the probe detuning. 
Note here that without a FORT the optimum cavity and 
atom detunings are not equal |p|~|ll|]. In our case, how- 
ever, the FORT effectively changes the atomic frequency 
in a position-dependent way and thus the precise value 
of the atomic detuning relative to the cavity detuning 
is largely irrelevant. Indeed optimum cooling conditions 
will exist in certain wells but not in others, which is one 
feature that allows one to distinguish various wells. 

Coupling the atom and the cavity to the remaining 
modes of the electro-magnetic field leads by a standard 
procedure to the master equation for the density operator 
of the coupled atom-cavity system, 



dp 
dt 



= —i[H, p\/fi ~ n{a + a, p) + Inapa^ 



{<T+a-,p} + — [ d 2 k V(d-e) 2 exp(-ifc-r) x 

87T J 

e 

X(7~p<r + exp(ifc ■ r), (5) 



with r the spontaneous decay rate and k the cavity decay 
rate. We are mainly interested in the strong-coupling 
regime, where ga 3> T,k. 

We treat the external (center-of-mass) degrees of free- 
dom of the atom c lassic ally, an approximation justified at 
the end of Section [V C . For a discussion of various inter- 



esting effects arising from the quantized external motion 
of an atom in a cavity QED field, we refer the reader to 
& 

In the quasiclassical approximation (i.e., where we re- 
tain the full quantum character of the internal degrees 
of freedom and of the cavity mode; see || for a full dis- 
cussion of this approximation), the integral in (^) can be 
evaluated to give the simpler result 



dp 
dt 



= —i[H, p]/h — n{a + a, p} + Inapa^ 



-{(7 + er ,p} + T<j pa + 



(6) 



The force acting on the atom consists of two parts, one 
due to spontaneous emission, whose mean vanishes on 
average, and the other part is represented by the operator 

F = -VH 

= -2hVS F {(J + (j- - 1/2) - hVg(a + a- + a+a), (7) 

which has contributions arising from the FORT poten- 
tial and from the interaction with the cavity mode. It 
was only the latter part that was considered in jlO| and 
that leads to 1-D cooling to temperatures of the order 
fcgT ~ min(7iK, KT/2). See also Refs. |l4|] for similar cal- 
culations on single atoms moving in cavity QED field, 
and Refs. jl5|,[l6) for calculations of diffusion of atoms in 
optical traps in free space. 

It can be shown [jlTj starting from a fully quantized de- 
scription, that the scmiclassical motion of the atom is de- 



scribed by a Fokker-Planck equation for the Wigner dis- 
tribution function containing (position-dependent) fric- 
tion and diffusion coefficients. Equivalently, we may use 
stochastic equations for the classical atomic position and 
velocity variables f and v of the form 



dr = vdt, 
(F) 

dv = —dt - (ivdt + BdW, 
m 



(8) 



where (.} denotes an expectation value, (3 is the friction 
tensor (with dimensions of a rate), m the mass of the 
atom, B is a tensor such that D = BB T jl is the velocity 
diffusion tensor (with dimension m 2 /s 3 ), and dW is a 
3-dimensional Wiener process which satisfies dWidWj = 
dijdt |Q. Starting with the expression (0) for the force 
operator, we can calculate (3 and D by the procedure 
outlined in the next Section. 



III. FRICTION AND DIFFUSION 

Refs. |To| employ Heisenberg equations of motion for 
various field and atomic operators to find friction and 
diffusion coefficients. These equations are not closed and, 
consequently, an approximation has to be made in order 
to find solutions. The natural assumption is to consider 
the weak driving limit (i.e., So <C k) and truncate the 
available Hilbert space to that part containing no more 
than a single cavity photon. This allows one to close the 
Heisenberg equations JTc| ] . Here we employ a different 
method (using the density matrix equations) to calculate 
friction and diffusion coefficients that does not require us 
to stay within the weak driving limit, but in addition we 
used Rcf. |l(J's procedure here to obtain results in the 
weak driving limit for verification purposes. In any case, 
it is still true that the most interesting regime is where 
only one or few photons are involved. Note that given 
the strong coupling between atom and cavity field, even 
a single photon is sufficient to lead to regimes far beyond 
the weak driving limit. In our examples we truncated the 
Hilbert space to photon numbers of around 4 or smaller. 
We refer to [[{J for an exposition on how to represent 
operators in truncated Hilbert spaces of precisely this 
form in a numerically convenient manner. 

The master equation (^|) is written as 



dp r 



(9) 



Numerically, the Liouvillian superoperator C is converted 
into a pre-multiplication operator by methods explained 
in |l9f| . In order to find friction and diffusion coefficients 
we apply a simple procedure, which yields these coef- 
ficients at zero velocity: this is sufficient for our pur- 
poses as the atom we are interested in, Cs, is relatively 
heavy. More precisely, the relevant dimensionless param- 
eters determining the velocity dependence of friction and 
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diffusion coefficients are kv/T and kvj n (see for instance 
pof ), and both are very small in all our simulations. In 
particular, T/k ~ 4.3 m/s and n/k ~ 3.4 m/s, while ve- 
locities in the trapping regime we are interested in (where 
atoms are localized in wells at low temperatures for times 
^> k~ 1 ,T~ r ) are around the Doppler limit velocity 



Vd 



hV/2 



8.8cm/s. 



(10) 



Also note that the standard procedure of continued frac- 
tions to calculate the full velocity dependence is not 
directly applicable to the present the poten- 

tial through which the atom is moving is not periodic 
(Af + A ). 

For an atom moving at velocity v we write 



d d r, 



(ii) 



and expand (g) in powers in v and solve for the steady 
state. The zeroth-order solution is then the steady state 
po at zero velocity: 

Cpv = 0, (12) 

while the first-order term p\ is determined by 

Cpi = v- V/9 - (13) 

The zeroth-order force is the steady-state force for an 
atom at rest, and is given by 



F Q = -Tr(poVtf). 



(14) 



Similarly, the friction coefficients follow from the first- 
order term in the force 



Pi = -Tr( Pl V#), 



by identifying 



(15) 



(16) 



where (3 is a 3-by-3 tensor. In our case ( 0), the gradi- 
ents along the cavity axis are larger in magnitude than 
those in the transverse directions by roughly a factor 
kwo w 150 (and around the cavity axis where the atoms 
spend most of their time the radial gradients are even 
smaller, of course). Since the friction coefficient scales 
with the product of two gradients (cf. Eqs (|l3|) and (pHl)), 
the largest element of the tensor f3 is the zz component. 
Next largest in magnitude are the off-diagonal compo- 
nents such as f3 xz and f3 zx . Their effects, however, can 
be safely neglected in our case: firstly, the force in the 
z direction proportional to —f3 zx v x is smaller than the 
friction force — (3 zz v z by roughly a factor kuiQ. Secondly, 
the force in the x direction —f3 xz v z is not a friction force 
(as it is not proportional to v x ) 1 and its contribution is 
averaged out because the oscillations in v z are faster than 



those in the x direction by another factor kwQ. Finally, 
the purely radial friction rates such as (3 XX are too small 
(<C Is -1 on average) to have any influence on the time 
scales considered here. Thus we take only (3 ZZ into ac- 
count. 

The diffusion coefficient, again at zero velocity, is cal- 
culated as follows. The standard method is to use the 
quantum regression theorem, and a particularly useful 
(for numerical purposes) interpretation of that theorem 
is given in jT~ 
given by 



The momentum diffusion tensor D p is 



D v — lim Re 

K t— >oo 



dr(F(t)F(t-r))-(F(t))(F(t-r)), 

(17) 



and its relation to the velocity diffusion tensor is D = 
Dp/m 2 . Before eliminating any degrees of freedom, the 
total system in fully quantized form is described by a 
time-independent Hamiltonian, which we denote by -fftot- 
In that case the time evolution of all operators is de- 
termined by exp(— iHtott), and two-time averages of the 
form (A(t)B(t — r)} as appearing in (17) can be written 
as 

(A(t)B(t - r)> = Tr [Aexp(~iH tot T)B Ptot (t) exp(ifl" tot T)] , 

(18) 

with ptot the density matrix of the total system. This 
expression formally contains the evolution of a density 
matrix over a time interval r starting from an initial 
density matrix p lni t = Bp to t{t)- The quantum regres- 
sion theorem now states that ([l8]) is still valid for the 
reduced density matrix that evolves under the Liouvil- 
lian C. That is, instead of (Eq) we may use 



(19) 



(A(t)B(t - r)) = Tr [Aexp(£ T )Bp(t)} 



In our case, £ is a time-independent operator and hence 
the right-hand side of (|l9|) can be evaluated by expanding 
exp(£r) in an exponential time series, as in the methods 
developed in . This then is the method we use here to 
evaluate the friction and diffusion tensors, and the results 
have been checked in the low-intensity limit by applying 
the different methods from |l(| to the same problem. 

Diffusion due to spontaneous emission is not obtained 
this way (as the bath of vacuum modes has been elimi- 
nated already), but can be obtained by standard meth- 
ods and gives an independent additional three compo- 
nents (-Dp)IP = Nih 2 k 2 r/2{<7 + a-) for i = x,y,z, with 
(.)o denoting a steady-state value and with the dimcn- 
sionless factor iVj depending on polarization. When the 
two-level system is formed by two Zeeman levels that are 
connected by circularly polarized light propagating in the 
z direction, we have N z = 2/5, and N x = N y — 3/10. 

Since the diffusion coefficients, just as the friction coef- 
ficients, scale as the square of a gradient, the largest com- 
ponent is D zz . Off-diagonal elements such as D xz and 
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D zx are, again, smaller by roughly a factor kw w 150, 
while the diagonal radial components such as D xx are 
in fact largely determined by spontaneous emission, and 
are of similar or larger magnitude than the off-diagonal 
elements. The proper way to take into account the off- 
diagonal elements of the diffusion tensor D is to diago- 
nalize D, and consider 3 independent diffusion processes 
along the axes of the basis that diagonalizes D with the 
eigenvalues of D as diffusion coefficients. Using the fact 
that D zz is large we can calculate both eigenvalues and 
eigenbasis perturbatively. The eigenvalues to first order 
are given by 



Dx'x' — D X x 

D z ' z ' = D zz 



D X 7.D 



D zz 
DxJh 
D„ 



where the . . . stands for terms of higher order in 1/ (kwo) 
while the axes change as 



z = z + x 



Dxz 

D zz 

. Dzx 



(21) 



The fact that z' is slightly tilted towards the x direction 
implies that a small part of the large diffusion coefficient 
D z i z i will contribute to diffusion in the x direction. This 
increase, however, is almost exactly compensated for by 
the decrease in D x / X i . In particular, the velocity in the x 
direction undergoes the following Wiener process: 



>2Dx 



2D 7 



D 2 
D 2 



AW. 



(22) 



In our case it turns out that D XX D ZZ ^> D zx (see Fig. ||), 
so that effects due to the off-diagonal elements of the 
diffusion tensor can in fact be neglected. The figure also 
shows that the previous considerations about the relative 
sizes of the various components of D do not just hold on 
average, but also locally. 




FIG. 3. For parameters to 

be used later, S /{2n) = 50MHz, and A p /(2tt) = -10MHz, 
iV e = 0.01 we plot here the ratio d = D 2 ZX /(D ZZ D XX ) as a 
function of position. 

Thus, friction is appreciable only along the cavity axis, 
while diffusion has two main contributions: from sponta- 
neous emission in all three directions, and a large diffu- 
sion along the cavity axis from fluctuations in the FORT 
and cavity QED forces. 



IV. NUMERICAL RESULTS 



The following results pertain to a Cs atom, with the 
ground state given by \&S\/ 2 ;F = 4;rriF = 4) and the 
(20) excited state by |6P 3 / 2 ;-F = 5;m,F = 5), so that Ao = 
852.4nm and r/(27r) = 5.2MHz. The cavity parameters 
are k/(2tt) = 4MHz and go/{2n) = 30MHz, and wq = 
20/i, which are typical for the experiments discussed in 
||. Furthermore, the values for Sq examined here are 
00/(2%) = 10, 50MHz. Both of these values are close 
to those explored in the actual experiment ||, and they 
contrast the behavior of atoms in shallow (So < go) and 
deep (So > go) wells. Typical values for N e range from 
10" 3 to 0.1. 



A. Dressed state structure 

We first focus on the atomic motion along the cavity 
axis. The simplest way to get a feeling for the results for 
(3 ZZ and D zz as a function of the probe detuning A p is to 
first consider the eigenenergies of the dressed atom-cavity 
states. When we neglect dissipation for the moment, and 
take the limit of no driving (N e — 0), we can easily find 
the energies of the lower dressed states \tp±) containing at 
most one excitation: the state containing no excitation is 
the ground state with an energy of Eq — — frSF(r), while 
the energies of the two dressed states in the manifold of 
states containing a single excitation are 



E± = huj a ± h^g{r) 2 + S F (r) 2 , 



(23) 



if the atom and cavity are on resonance. The excited 
dressed states are given by 



IV'-) = sin%, 1) + cos0|e,O), 



(24) 



with 



g 2 + (vV + sJ - s F y 



cos 8 



9 2 + (V^+s 2 - S F ) 2 



(25) 



In Figures g (10MHz FORT) and | (50MHz FORT) we 
plot the transition frequencies (relative to uj a ) from the 
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ground state to these two excited states as functions of 
position, i.e., 

A± = S F (r) ± vW + W' (26) 

This expression along with the figures explicitly shows 
that the main features of the atom-cavity system are de- 
termined by the ratio So/ go. It furthermore shows an 
important difference with the situation of trapping with 
a FORT in free space. The fact that the excited state 
shifts up while the ground state shifts down implies that 
ground and excited states are trapped in different posi- 
tions in free space. In the presence of the quantized cavity 
field, however, both the lower excited dressed state and 
the ground state are now shifted down. This may improve 
trapping and cooling conditions, as detailed below. 
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position [X] 

FIG. 4. Transition frequencies A± relative to the bare 
atomic frequency from the ground state to the lower two ex- 
cited dressed states as functions of the position of the atom 
along the cavity axis (i.e. p = 0) . Here So/(2n) = 10MHz. 
Also indicated by the dotted line is the probe detuning used 
in Fig. | A p /(2n) = -28MHz. 
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FIG. 5. Same as previous figure but for So/(2tt) — 50MHz. 
Also indicated by the dotted line is the probe detuning used 
in Figure ^, A p /(2tt) = -10MHz. 

B. Cooling mechanisms 

We now take a closer look at cooling mechanisms. In 
the regime of weak driving, we will find that the friction 
coefficient j3 zz is positive (corresponding to cooling) when 
the probe field is tuned slightly (by an amount ~ n, T/2) 
below the transition to the relevant dressed state while 
for blue detuning the friction coefficient is negative, lead- 
ing to exponential heating of the atom's velocity. This 
can be understood by analogy with Doppler cooling: by 
tuning below resonance, the process of stimulated ab- 
sorption followed by spontaneous emission leads to a loss 
of energy, while the maximum cooling rate is achieved 
by maximizing the product of excitation rate and detun- 
ing. Now looking back to Figs. |] and || one sees that 
the variation of A+ with position is larger than that of 
A_, because both the ground state and the lower excited 
dressed state shift down, while the upper excited dressed 
state shifts upward. Generally speaking, for cooling pur- 
poses it is better to tune to the lower excited state so as to 
have smaller spatial variations in cooling rates. More im- 
portantly, the upper excited state energies decrease with 
increasing radial distance, whereas the lower excited state 
energy increases. Thus, for the upper state the probe de- 
tuning changes from red to blue, so that an atom cooled 
on axis will in fact be heated if it moves away radially. 
For the lower dressed state the probe detuning becomes 
more red, so that an atom that is optimally cooled on axis 
will still be cooled away from the axis, but at a lower rate. 

The most popular explanation for intra-cavity cool- 
ing (!(]] exploits analogies with Sisyphus cooling 
although another explanation for cavity-based cooling 
based on asymmetries in coherent scattering was recently 
put forward in |y]]. Here we illustrate the Sisyphus pic- 
ture for cooling inside optical wells within an optical res- 
onator, using a very simple dressed-state picture, that 
makes use of only the lower dressed state and the ground 
state, relevant in the low-excitation limit. We choose 
one particular well, from z = 2.0\p to z = 2.5Xf, and 
one particular set of parameters given in the caption of 
Fig. ||. In that figure we plot the decay rate 7_ of the 
lower dressed state and the excitation rate from ground 
to the dressed state, as functions of position. In the 
weak driving limit the decay rate is given by 

7_ = (^_|Ka + a + rcr + cr _ /2|^_) 

= sin 2 6k + cos 2 6T/2, (27) 

and the excitation rate by 

0_ = |(<7,0|£ (a + + a)h/'-}l 

= £ o |sin0|. (28) 
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These two quantities, together with the detuning of the 
probe field from (dressed-state) resonance determine the 
steady-state population in the lower dressed state, ac- 
cording to 



n 2 



(A_-A p )2+ 7 2 



2 ■ 



(29) 




2.2 2.3 
position [X F ] 

FIG. 6. Decay rate and excitation rate of the lower dressed 
state as functions ofposition along the cavity axis (p = 0). 
Here and in Figures [?j and ^| we chose the following parame- 
ters: N e = 0.001, So/(2n) = 50MHz, A p /(2tt) = -10MHz. 

The population n_ is plotted in Fig. |t], along with 
the transition frequency A_. These two quantities are 
sufficient to understand the Sisyphus cooling mechanism. 




2.2 2.3 
position [X ] 

FIG. 7. Transition frequency to and populations in the 
lower dressed state as functions of position along the cav- 
ity axis. Note that the equilibrium position of the atom is 
around z = 2.25Af- 

Since an atom in the ground state is moving in a 
conservative potential around the equilibrium position 
z = 2.25Af, the following Sisyphus picture should be 
taken as to apply to the motion of the atom in addi- 



tion to that conservative motion (see (p0[)). Suppose, for 
example, that the atom is at position z = 2.2\p and 
moving towards the right (cf. Fig. ^) . The probability to 
be in the excited state now decreases (according to the 
lower part of Fig. Q), while the energy of the excited state 
relative to the ground state is increasing: in other words, 
an atom in the excited state is climbing uphill (again, 
in relation to the ground state), but will likely make the 
down transition to the ground state, thus leading to cool- 
ing at that particular position. Similarly, at z = 2AXp 
an atom moving to the left is going uphill while having an 
increased chance of decaying to the ground state, again 
leading to cooling. This picture in fact shows that the 
cooling rate is expected to be proportional to the gradi- 
ent of n- and the gradient of A_. More precisely, the 
force on the atom at position z is approximately given 
by 



. (ISf 
dz 

dS F 
dz 



Tin_(z - 
hri-(z) 



- v/7- 
dA_ 
dz 



dz 

hv dri- dA_ 



7- 



dz dz 



(30) 



where the argument of n_ indicates the lag between the 
atom reaching a position z and reaching its steady state, 
with the lag time scale determined by the inverse decay 
rate from the dressed state. From the second line we see 
that the friction coefficient f3 zz is approximated by 



R : 



h dri- dA_ 



7717- 



dz dz 



(31) 



Indeed, Fig. || shows the similar behavior of /3 ZZ and R 
as functions of position. 
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FIG. 8. Cooling rate /3 ZZ (solid curve) and the product of 
gradients of dressed-state population and transition frequency 
R (dashed curve) as functions of position along the cavity axis. 
The similarity between the two curves confirms the validity 
of the Sisyphus cooling picture. 
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C. Friction, diffusion and equilibrium rms velocities 

In Figures p|-|lO| we give examples of friction and dif- 
fusion coefficients for both the 10 and 50 MHz FORTs, 
as functions of the atomic position. They illustrate the 
point that in the low-excitation limit red (blue) detun- 
ing leads to cooling (heating) (cf. Figs. [| and g). They 
moreover clearly show how all wells are quantitatively 
different, with cooling rates and diffusion strengths dif- 
fering by orders of magnitude over the various wells, and 
with (3 ZZ being negative in some wells, and always positive 
in others. This of course also implies that the tempera- 
tures reached by atoms in thermal equilibrium vary with 
position. 
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FIG. 9. Friction and velocity diffusion coefficients /3 ZZ and 
D zz as functions of the atomic position (in units of Af) along 
the cavity axis. Here N e = 0.001, So/(2ir) = 10MHz, and 
A p = -28 x 2ttMHz. Cf. Fig. f§ 

For the case of the shallow FORT we consider weak 
driving (N e — 0.001), whereas for the deeper FORT the 
driving field is taken to be stronger by an order of mag- 
nitude. The stronger driving field increases cooling rates 
while the fact that deeper wells trap the atoms better 
means that correspondingly larger diffusion rates still can 
be tolerated. 
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FIG. 10. Same as previous Figure but for 



4 

0.01, 



5 /(2tt) = 50MHz, and A P 



-10 x 2ttMHz. 



The stable equilibrium points z^ are located around 
the maxima of Sp, i.e. around z n = (n — 1/2) Air/2 for 
integer n, because it is the FORT that gives the main 
contribution to the total force (even for the smallest value 
of So = 2tt x 10MHz considered here). The cavity QED 
field gives only a small correction to the force and hence 
to the equilibrium position. In each equilibrium point, 
we can define a measure for the expected rms velocity 
of the atom along the z axis in thermal equilibrium by 
considering averages over local wells 



A. 



if A, >0, 



(32) 



in terms of the friction and diffusion coefficients. This 
averaging procedure gives a sensible measure for the rms 
velocity only if the atom indeed samples the whole well. 
This condition is fulfilled for the relatively shallow wells 
originating from So = 2tt x 10MHz, and Fig. [ll] uses this 
averaging procedure. For the 50MHz FORT, however, we 
averaged over only part of the well, namely a region of 
size Af/10 symmetrically around the equilibrium point. 
This choice is rather arbitrary, and thus Fig. [l2]just gives 
an indication of what rms velocities to expect for atoms 
trapped in the corresponding wells, although the simula- 
tions in fact do confirm these values. 
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FIG. 11. The values of Vr ma in the eight equilibrium points 
as defined in (|3^) by averaging over the entire well. In all 
cases So = 2-zr x 10MHz. Triangles correspond to a detuning 
A p /(2tt) = -28MHz, squares to A p /(2tt) = -23MHz, and 
circles to A p /(2n) = 15MHz . Note the points on the latter 
curve on the axis indicate that the friction coefficient is nega- 
tive, so that there is in fact no cooling and v* ms is not defined. 
They do not indicate cooling to u rms = 0. 



peratures achievable are determined by the Doppler ve- 
locity Vd- More precisely, the lower limit on rms veloci- 
ties along the cavity axis is expected to be 

V D = f^, 03) 

where the factor 0.7 = (1 + 2/5)/2 comes from the fact 
that in our case the diffusion due to spontaneous emission 
in the z direction is 2/5th of the full 3-D value. Wc 
tested that for smaller k the rms velocitie s indee d do 
become even smaller, now determined by ^Tin/m, thus 
confirming predictions of |Tc| ] . 

Finally, we note that the quasiclassical approximation 
used throughout this paper is justified as neither the re- 
coil limit is reached nor the resolved-sideband limit, i.e. 

TiT/2 > (hk) 2 /m, 

fiT/2 > hu osc , (34) 

with v oac the oscillation frequency of the atom in a well 
(see below), although in some cases the latter condition 
is only marginally fulfilled, namely when v osc = 600kHz, 
which is only a factor 4 smaller than r/(47r). 
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FIG. 12. As previous Figure, but for So = 2ir x 50MHz and 
N e = 0.01, and where the average is taken over a region of size 
Af/10 around the equilibrium point. The probe detunings 
were A p /(27r) = —10, —5, 100MHz for the squares, triangles, 
and circles, respectively. Note the points on the axis indicate 
that the friction coefficient is negative, so that there is in 
fact no cooling and ii, 2 ,,, is not defined. They do not indicate 
cooling to w rms = 0. 

We see here that depending on the probe detuning, 
the atom will be cooled to low temperatures either in all 
wells, or only in wells where g is large in the equilibrium 
point, or only in wells where g is small. This shows the 
flexibility that a FORT beam adds: one can predetermine 
to a certain degree in which well the atom will be trapped 
(and cooled) for longer times and in which it will not be. 

Under the current conditions n > T/2 the lowest tem- 



D. Saturation behavior 

We now briefly turn to the question of the nonlinear 
behavior of the atom-cavity system with increasing exci- 
tation. In the absence of saturation effects, both friction 
and diffusion coefficients would increase linearly with N e . 
For the same parameters as Fig. |9[ Figure |lj shows non- 
linearities setting in around N e — 0.01. The friction coef- 
ficient even starts to decrease around N e = 0.1 as a result 
of the local values of f3 zz becoming negative where they 
were positive in the weak driving limit. The concomitant 
effect on the v* ms is shown as well. 
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FIG. 13. The average values of D zz and /3 ZZ as functions 
of the driving field strength N e in the well extending from 
z = 2Af to z = 2.5Af, for the 10 MHz FORT, where 
A p = —28 x 27rMHz. In the lower part the corresponding 
values for the rms velocity v^ ms are plotted as a function of 



E. Simulations 

We also performed Monte-Carlo simulations of the 3-D 
motion in given wells by solving the Langevin equations 
(0) for position and velocity (see also The experi- 

mental procedure switches the FORT field on only when 
an atom has been detected and when it consequently has 
partly fallen through the cavity already We accord- 
ingly fix initial conditions as follows: We start the atom 
on the cavity axis, and we fix the downward velocity to 
be v x = 10 cm/s. Furthermore, we chose v z = cm/s, 
and the initial position along the z axis to be Xp /8 away 
from the equilibrium point. The initial position and ve- 
locity were fixed so that all variations in trapping times 
and rms velocities are solely due to the random fluctua- 
tions of the forces acting on the atom, rather than from 
random initial conditions. Experimentally these two are 
mixed of course. 

Since atoms with these initial conditions do not possess 
angular momentum around the z axis, this in some sense 
represents a favorable case (although the atoms are not 
put in the bottom of the well). However, in the course of 
their evolution the atoms do acquire angular momentum 
so that this is in fact not a severe restriction. For more 
detail see below (Figure ||^) . 

In Figure [l4]we plot the results of simulations of 1000 
trajectories for an atom in the shallow well of 10 MHz. 
We plot the average rms velocity along the cavity axis 
as a function of trapping time for each trajectory. Here 
we defined the "trapping time" as the time spent by the 
atom in one particular given well of size A^/2. The actual 
trapping time inside the cavity may be longer, obviously, 
as the atom may subsequently get trapped in different 
wells. For very short trapping times, v rms is determined 
by the initial condition, but for longer times lower tem- 
peratures corresponding to those calculated in Fig.|ll| are 
reached. Note however that the simulations were done in 
3-D, and as such do not necessarily give the same temper- 
atures as predicted for on-axis (1-D) motion in Figures [TT| 
and |lj. Nevertheless, the effect of the atoms' radial mo- 
tion is apparently not strong, and in fact atoms leave the 
well while still being trapped radially. This is partly due 
to the fact that all (especially heating) rates in the radial 
direction are smaller by a factor kwo ~ 150 than those 
in the axial direction. 
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FIG. 14. In the upper part of the plot each data point 
gives the trapping time and the average u rms resulting from 
a single trajectory. Identical initial conditions were chosen 
for each trajectory: each atom started at z — 2.125Af with 
v z = and v x = lOcm/s (downward). Other parameters 
were: iV e = 0.001, A p = -28 x 2ttMHz, So = 10 x 27rMHz. 
The lower part gives a histogram of the probability P(T) for 
an atom to be trapped longer than a time T. A fit of the 
tail of this distribution to an exponential oc exp(— T/r) gives 
r = (25 ± 2) ms 

About half of the atoms is basically not trapped at 
all. The remaining atoms have a probability P(T) to be 
trapped longer than a time T, with P(T) decaying ex- 
ponentially with T. The average trapping time for these 
parameters is found to be r ~ 25ms, as shown in Fig. 

In Figures |lj and [l6] we plot for the same 10MHz 
FORT an example of a single trajectory, after the atom 
has spent 4ms in the trap. The oscillation frequencies 
along the z and the radial directions differ by two or- 
ders of magnitude (since kwo m 150): in the z direc- 
tion the oscillation rate is ~ 200kHz, in the radial di- 
rection ~ 2.2kHz. The photon transmission follows both 
these oscillations so that in principle the atomic motion 
in both axial and radial direction is detectable. Exper- 
imentally, though, the oscillations along the cavity axis 
may be too fast to be accessible. In particular, the aver- 
age rate at which photons leaking out through one end 
of the cavity are detected is at most (the efficiency is 
less than 100%) equal to the cavity decay rate multi- 
plied by the average number of photons inside the cavity. 
For the parameters of Fig. [l^ this amounts to a rate 
~ 0.01 x k ~ 2.5 x 10 5 /sec, which corresponds to just 
about one photon per oscillation period. 

The figures show that when the atom is in a position 
where it is not coupled to the cavity {g = 0) , the number 
of photons in the cavity drops to N e — 0.001. Similarly, 
when the atom moves away radially, the transmission 
drops. 
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FIG. 15. Snapshot of a single trajectory, with parameters 
as in Figure [L4|. The upper plot gives the z coordinate of the 
atom as a function of time, the lower plot gives the transmis- 
sion (in fact the number of photons (a + a) inside the cavity) 
in that same time interval. Note the time scales here differ 
by two orders of magnitude from those of Figure um 
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FIG. 16. For the same trajectory as the previous figure, the 
upper plot gives the radial distance to the cavity axis, p in 
units of Wo as a function of time, the lower plot gives again the 
number of photons inside the cavity during that same time 
interval. The atom has a nonzero angular momentum along 
z and does not cross the z axis. 

To make a direct comparison with the trapping times 
achieved in the experiment we now turn to the case 
of a 50 MHz FORT. We plot rms velocities vs trapping 
times for 300 trajectories for an atom trapped in the well 
ranging from z — 2\p to z = 2.5\p. 
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FIG. 17. As Fig. |ll| but for N e = 0.01, A p = -10x2ttMHz, 
So — 50 x 27rMHz. The mean trapping time is r = 250±20ms. 

For the parameters of Fig. [l?] the atom is either 
trapped for long times (> 10ms) or only for a short time 
(< lms), both with about 50% probability. In the latter 
case the rms velocity is determined just by the (arbitrar- 
ily chosen) initial condition and is around 30cm/s, but 
for longer trapping times the effects of cooling are visi- 
ble. Thermal equilibrium is reached with v rms ~ 8cm/s, 
thus confirming the results of Fig. [l2]. The distribution 
of trapping times again follows an exponential law, and 
the average trapping time, as determined from the tail of 
the distribution, is r w 250ms, which is ten times longer 
than for the (fluctuating) 50MHz FORT used in @. This 
shows the great potential of holding single atoms in the 
cavity for extended periods of time if the intensity fluc- 
tuations of the FORT beam can be minimized. Experi- 
mental efforts along this path are currently underway. 

Also for this case we plot snapshots for a single trajec- 
tory, taken after the atom has spent 25ms in the trap. 
Compared to the 10MHz FORT, the oscillations of the 
atom along the cavity axis and in the radial direction be- 
come faster by about a factor of 3. The axial oscillation 
frequency is about 600kHz, while along the radial direc- 
tion the oscillations occur at a rate 6.2kHz, i.e., again 
slower by two orders of magnitude. In this case, the 
photon transmission still follows directly the axial os- 
cillations but no longer follows the radial excursions of 
the atom, as now the fluctuations in the magnitude of g 
at the atom's position along the cavity axis are in fact 
larger than those due to the radial excursions of the atom. 
This is partly due to the fact that in the simulations here 
the driving field is stronger than for the 10MHz example 
above so that fluctuations in the atomic motion occur at 
a shorter time scale, and partly simply because the radial 
excursions are small. Fig. |l9|(c) shows that it is primar- 
ily the axial fluctuations that determine the variations in 
the numbers of photons inside the cavity. 

Generally speaking, the axial excursions determine (lo- 
cal) minimum and maximum transmission levels (as in 
Fig. 15). When these minima and/or maxima depend 
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on the radial position, then the radial motion could in 
principle be visible in the cavity transmission level. This 
depends in turn on whether the axial fluctuations on the 
time scale of the transverse motion are sufficiently small 
so as not to hide the radial dependence. There seems 
to be no simple general rule how this interplay between 
radial and axial motions depends on detunings, driving 
strength, and the particular well. 
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FIG. 18. Snapshot of a single trajectory, with parameters 
as in Figure The upper plot gives the z coordinate of 
the atom as a function of time (with the atom released with 
standard initial conditions at t — 0) . The lower plot gives the 
transmission (in fact the number of photons inside the cavity) 
in that same time interval. Note the time scales here differ 
by two orders of magnitude from those of Figure |l9[ 



1.0.05- 




25 



25.2 



25.4 



25.6 



25.8 



26 



0.03 



6.02 - 





25.4 25.6 
time [ms] 

FIG. 19. For the same trajectory as the previous figure, (a) 
the radial distance to the cavity axis, p in units of wq as a 
function of time, (b) the number of photons inside the cavity 
during that same time interval, and (c) the position along the 
cavity axis in units of \f- 

In contrast, in a different well, the one ranging from 
z = 0.0 to z = 0.5Af, the photon number in the cavity 
does follow the radial motion, as the radial excursions 
become larger. Perhaps more importantly, the average 



transmission level is higher by more than a factor 2 com- 
pared to the previous case, as a result of g being larger 
in this well (cf. Fig. |^). This shows how, in principle, 
different wells may be experimentally distinguished via 
the transmission of the probe field through the cavity. 
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20. For an atom in the well ranging from z — 

-10 x 2ttMHz, and 



A F /2, for N e = 0.01, A p 



to z 

So — 50 x 27rMHz, the upper plot gives the radial distance to 
the cavity axis, p in units of wo as a function of time, the lower 
plot gives the number of photons inside the cavity during that 
same time interval. 

We also simulated the motion of an atom trapped un- 
der more adverse conditions, namely for an atom in the 
well [z = Aj? — > 1.5Af] at a probe detuning A p /(27r) = 
— 5MHz. According to Fig. [l^, the atom is not cooled 
on axis under these conditions (i.e. the average friction 
coefficient around the equilibrium point on the z axis is 
negative). This is confirmed by Fig. the mean trap- 
ping time for an atom starting at z = 1.125Af is now 
very short, about 1.6ms, while the average rms velocity 
is uf 



condition. 



28 cm/s, as determined essentially by the initial 




0.002 



0.004 0.006 
time [s] 



0.008 



12 



FIG. 21. As Fig. but for N e = 0.01, A p = -5x 2ttMHz, 
So = 50 x 27rMHz. The initial position of the atom is 
z — 1.125Af. The mean trapping time is r = 1.6 ± 0.1ms. 

Finally, we consider the influence of different initial 
conditions on trapping and cooling. All the results so 
far were obtained by considering atoms that initially are 
moving on axis. Thus, they have no angular momen- 
tum along the z axis, nor any radial potential energy. 
Figure shows a plot of rms velocities vs. trapping 
times for atoms trapped under the same conditions as 
for Fig. [I?]) (i.e., in the well from z = 2A to z = 2.5A, for 
A p = -10 x 2ttMHz, N e = 0.01, and S = 50 x 2ttMHz), 
but with different (nonzero) values for the initial angular 
momentum. 
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FIG. 22. Rms velocities vs. trapping times for atoms 
trapped under the same conditions as for Fig. [l?] but with 
different initial radial conditions for y. In particular, for plot 
(a) the initial conditions on y is y = 0.2wn. for (b) y = 0.5wn 
and for (c) y = idq. Since v x = — lOcm/s, the atoms have dif- 
ferent angualr momenta along z in these cases, and different 
initial potential energies. 



F. A different trapping structure 

We now consider a different case where the atomic ex- 
cited state is assumed to be shifted down by the FORT 
field, just as the ground state is (see, for instance jf2|). 
This can be achieved by using a FORT that is (red) de- 
tuned in such a way that the excited atomic state is rel- 
atively closer to resonance with a higher-lying excited 
state than with the ground state. This situation at first 
sight looks even more appealing for trapping purposes, 
as now both excited and ground state will be trapped in 
the same positions. Moreover, fluctuations in the force 
due to the FORT are dimished. 

We consider only the 50MHz FORT here, and compare 
this case to the previous 50MHz FORT case, and in par- 
ticular we refer the reader back to Figs. [|, [jl], and |T^. 
For ease of comparison we keep A^ the same, and assume 
for simplicity that the excited state is shifted down by an 
amount Sf, so that the shifts of the ground and excited 
state are in fact identical. 

The fact that ground and excited state see the same 
potential, implies that the transition frequencies to the 
dressed states are simply periodic in space with period 
Aq, as shown in Fig. p3l rather than aperiodic as in Fig. |^. 




Obviously, the more initial potential energy the atom 
has, the less likely it is to be trapped. In fact, the angular 
momentum does not play any role here, as confirmed by 
similar calculations with initial conditions chosen such 
that the atoms have no initial angular momentum but 
have the same potential energy. The results are the same 
in that case. For atoms starting at y — 0.2u>o the trap- 
ping times and rms velocities are basically not affected, 
and the trapping time is still around 250ms. But for 
atoms starting at y = 0.5i«o the effect of their increased 
potential energy leads to clearly shorter trapping times 
(by roughly a factor of 2), and for atoms starting at 
y = Wq this effect is even more pronounced with a de- 
crease in trapping time of about a factor of 10. 
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FIG. 23. Transition frequencies A± for the case where the 
atomic excited state is assumed to be shifted down by the 
FORT field by the same amount Sf as is the ground state, 
for the 50MHz FORT. A detuning of -35MHz is indicated by 
the dotted line. 

Similarly, the fluctuations in the force due to the FORT 
now vanish, as both ground and excited state undergo 
the same shift, so that the diffusion coefficient is periodic 
with period Ao- Also the friction force arises only from 
the cavity QED part and is periodic. Yet, the different 
wells are not equivalent. The forces are, just as before, 
driven by both cavity QED field and the FORT, and the 
value of g at the antinode of the FORT still varies over the 
different wells. This is illustrated in Fig. pi| where the rms 
velocities in the 8 different wells are shown, along with 
friction and diffusion coefficients. Since in this example 
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the probe field is detuned below the lower dressed state, 
one has cooling everywhere in space. 
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FIG. 24. Friction and diffusion coefficients, and the re- 
sulting rms velocity as functions of position al ong t he cav- 
ity axis for the trapping structure of Section [VF. Here 
A p /(2tt) = 



-35MHz, and iV e = 0.01. 



The simulations show that the mean trapping time is 
smaller, although the rms velocities are just as small as 
before. The reason is the less favorable cooling condition 
away from the cavity axis. In particular, for the param- 
eters used here the expected rms velocity v* ms steadily 
increases to 90cm/s at a radial distance p = 2wq, while 
for the simulations of Figs. [Tt], v* ms is increasing only 
slowly to 12cm/s. This large difference can be under- 
stood by noting the difference in dressed state structures 
between the two cases. For the case of Fig. [| the tran- 
sition frequency to the lower dressed state around the 
equilibrium position z « 2.25\p does not change much 
with increasing radial distance, so that the probe field in 
that trapping region is always detuned below resonance 
by an amount that stays more or less constant. For the 
dressed state structure of Fig. |23|, however, the probe 
detuning increases from > 5MHz to > 35MHz below res- 
onance, thus leading to much worse cooling conditions. 
In other words, the presence of opposite level shifts due to 
the FORT makes the spatial variation of the transition 
frequency to the lower dressed state smaller, compare 
A_ = Sf — V 'S'p + g 2 to A_ = — g, especially when 
g < S F - 

The alternative trapping potential is, therefore, not 
necessarily more favorable for trapping purposes. On the 
other hand, all atoms are captured now and are trapped 
for at least 10ms. This can be understood from the simple 
fact that here the friction coefficient is positive in the 
entire well. 
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FIG. 25. As Figur e but for the different trapping struc- 
ture of Section [V Fj . The initial position was z — 0.125Af, 
and further parameters were N e = 0.01, A p = —35 x 27rMHz, 
So — 50 x 27rMHz. The mean trapping time is r = 28 ± 2ms. 



V. SUMMARY 

We analyzed cooling limits and trapping mechanisms 
for atoms trapped in optical traps inside optical cavities. 
The main distinguishing feature from previous discus- 
sions on cooling of atoms inside cavities is the presence 
of the external trapping potential with a different spatial 
periodicity as compared to the cavity QED field. This 
not only provides better cooling and trapping conditions 
but the different spatial period makes the various poten- 
tial wells qualitatively different. Atoms can be trapped 
in regions of space where the coupling to the cavity QED 
field is maximum, minimum or somewhere in between. 
Depending on the laser detuning, cooling may take place 
only in wells where the atom is minimally coupled to 
the cavity QED field, or where it is maximally coupled. 
This allows one in principle to distinguish to a certain de- 
gree the different atomic positions along the cavity axis, 
namely, by comparing 

• the average transmission level 

• the fluctuations of the cavity transmission 

• the total trapping time 

which reflect, respectively, the average atom-cavity cou- 
pling, the temperature of the atom and under certain 
conditions the radial motion, and the overall cooling and 
trapping conditions. This is an important additional tool 
useful for eventual control of coherent evolution of the 
atomic center-of-mass degrees of freedom, as relevant to 
performing quantum logic operations. 
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